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We show formally how the application of boundary layer techniques to 
a well-stirred, concentrated, multicomponent electrolyte yields essentially 
the same results as those for a dilute, binary electrolyte. In particular, 
the coupled integral equations for the reactant concentration and the 
electric potential, derived by J . L. Blue for a dilute, binary electrolyte, 
remain valid in the general case with modified values of certain combina- 
tions of physical constants appearing in the equations. 

I. INTRODUCTION 

In this paper, we use the method of matched asymptotic expansions 
to derive boundary layer equations for a well-stirred, concentrated, 
multicomponent electrolyte containing a single ion reacting at the 
electrodes. We have two main objectives. First, we wish to show how 
the boundary layer approximation, up to now applied mainly to dilute, 
binary electrolytes or to the case of excess supporting electrolyte, 1 
may be systematically generalized. Second, and more important, we 
show how the generalization leads to equations almost identical with 
those valid for a dilute, binary electrolyte, for which J. L. Blue 2 has 
derived an elegant and efficient method of solution. 

In the general case, the singular perturbation in the boundary layer 
leads to a set of coupled convective diffusion (cd) equations for the 
concentrations or, in other words, a vector CD-equation in place of the 
single, scalar CD-equation for a dilute, binary electrolyte. The vector 
CD-equation is diagonalized by introducing the eigenvalues and eigen- 
vectors of a reduced diffusivity matrix. The resulting uncoupled, scalar 
CD-equations can then be solved by Blue's method. 

The perturbation yields linear CD-equations in the electrolyte, but 
the strongly nonlinear dependence of the electrode current density on 
the electrolyte potential and the reactant ion concentration is retained. 
In fact, the first approximation is simply the so-called "secondary 
current distribution," the potential satisfying Laplace's equation in 
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the electrolyte, but satisfying nonlinear boundary conditions at the 
electrodes. Furthermore, if, as is often the case, the current-carrying 
capacity of the electrode surface reaction is large compared with that 
due to ohmic conduction, the potential is approximately constant on 
the electrodes, giving the so-called "primary current distribution." 

We begin in Sections II and III with a semiqualitative discussion 
of some of the approximations involved. Sections IV through IX 
contain the detailed analysis. The reader not interested in the details 
may skip to Sections X and XI, which contain very brief summaries 
of the formal perturbation solution and the equivalent boundary 
integral solution developed by J. L. Blue. 

II. THE FIRST APPROXIMATION 

A first approximation to the cathode current density serves as a 
convenient introduction to the general concentrated, n-component 
electrolyte. In a well-stirred electroplating bath, the ion concentrations 
are almost constant, except in thin diffusion boundary layers surround- 
ing the electrodes. If the electrode current densities are bounded, so 
are the concentration gradients in the boundary layers. Thus, the 
changes in concentrations across the boundary layers are small when 
the boundary layer thicknesses are small enough. To a first approxi- 
mation then, the concentrations are constant throughout the elec- 
trolyte, current being carried entirely by ohmic conduction. In the 
electrolyte, the dimensionless potential ip (equal to FQ/RT, for Fara- 
day constant F, gas constant R, and absolute temperature T) satisfies 
Laplace's equation 

VV = 0. 

On the electrodes, the current densities are determined by surface 
reactions producing cations (e.g., CU++) at the anode and depositing 
them at the cathode, the reaction rates being functions of the local 
concentrations and potentials. The potential <p then satisfies the 
boundary condition 

{— j a , on the anode, 
j e , on the cathode, 
0, otherwise, 

where n is the outward normal, fc (0) is the dimensionless bulk con- 
ductivity, and, for e = a, c, 

je = j* (0) [exp(a e | <p - <P„\) ~ exp(-/3 e | ip - <p e \)J 

The dimensionless ion exchange current density j'* (0) , here constant, 
but in general a function of the concentrations, and the exponents a e , 
/3« characterize the surface reaction rates. The quantity <p — <?« is the 
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so-called "electrode surface overpotential" for electrolyte potential <p 
and electrode potential <p e - 

When j* l0) is large, | <p — ip e \ must be small or, to a first approxi- 
mation, <p = ip e on the electrodes, giving the so-called "primary current 
distribution." Note, however, that the primary current distribution 
leads to unbounded current density at electrode edges or reentrant 
corners, while the secondary current distribution does not. The ex- 
ponential form of the electrode current relation, however, suggests 
that the primary and secondary current distributions differ only in a 
very small neighborhood of such a boundary singularity. 

Once <p has been determined, the electrode current densities can be 
calculated. Furthermore, the approximation does not depend on the 
specific electrolyte, but only on the existence of thin boundary layers 
across which the concentration changes are small. It thus applies to 
dilute, binary and concentrated, multicomponent electrolytes alike, 
although, of course, the various constants, e.g., the conductivities, may 
be vastly different in the two cases. 

To obtain a better approximation, we must include the effect of 
concentration variation across the boundary layers. Before doing so, 
however, let us consider a method of practical solution for <p, the 
method of boundary integrals widely used for the numerical solution 
of potential problems. Since <p is harmonic, it satisfies Green's third 
identity, which in the plane has the form 

2*^ = / p [5j (WW q) - *w>H <*» «)] *<g>» 

where r is the curve bounding the electrolyte and G(P, Q) 
~ In (1/| PQ | ), as the point P — > Q, and is harmonic otherwise. When 
P lies on r, the identity becomes an integral equation for (p, nonlinear 
because d<p/dn is a nonlinear function of <p. Very efficient numerical 
methods for solving this equation have been developed. 

III. THE BINARY ELECTROLYTE 

For the dilute, binary electrolyte, J. L. Blue 2 has extended the above 
boundary integral method to include the effect of concentration changes 
across the boundary layers. He derived a relation giving the concen- 
tration changes across the boundary layers in terms of integrals of 
the product of electrode current density and a kernel function whose 
form depends upon the local velocity distribution. These concentra- 
tion changes affect the electrode current density in two ways. First, 
as we have already noted, the ion exchange current density is a func- 
tion of concentration. Second, the potentials at the electrodes differ 
from those at the boundary layer edge, where the previous <p is actually 
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evaluated, by the so-called "concentration overpotentials" A<p e , 
e = a, c. These overpotentials are due to the net current carried by 
ions diffusing at different rates through the boundary layers. To the 
present degree of approximation, they are linear functions of the 
concentration changes across the boundary layers. 

In a binary electrolyte, in which the two ion concentrations must be 
proportional, since the electrolyte is electrically neutral, A<p e is pro- 
portional to the integral giving the concentration change. This con- 
stant of proportionality, the constant giving the concentration to be 
inserted into the ion exchange current density, and the bulk conduc- 
tivity are the only constants containing material properties in the final 
relations. They thus characterize the electrolyte completely in this 
case and may thus be used in place of the individual diffusivities and 
mobilities. 

In the above, we did not qualify the term "binary electrolyte" with 
the term "dilute." The discussion of the previous paragraph applies 
as well to concentrated binary electrolytes, provided that the boundary 
layer concentration changes are small. In this case, in the boundary 
layers, the equations are linearized around a state of constant bulk 
concentrations, so that the diffusivities, in general forming a 2 X 2 
matrix, and the mobilities, both functions of the concentrations, are 
constants, evaluated at the constant values of the bulk concentrations. 
Using the charge neutrality condition, the potential and one concen- 
tration are eliminated from the two boundary layer equations govern- 
ing the ion fluxes. The resulting equation has exactly the same form 
as that for a dilute, binary electrolyte, so that the solution has the 
same form, except that the proportionality constants have different 
values. 

In the following, we extend this result to concentrated, multicom- 
ponent electrolytes, the case of the ternary (w = 3) electrolyte being 
of particular interest. As above, we linearize the equations around a 
state of constant bulk concentrations, using the boundary layer 
approximation (not the assumption of excess supporting electrolyte). 
We eliminate the potential and one concentration to obtain a set of 
(n — 1) boundary layer equations for the concentrations, coupled by 
a reduced (n - 1) X (n — 1) diffusivity matrix. This system is solved 
by diagonalization, using the principal values and vectors of the re- 
duced matrix. Finally, these solutions are combined to calculate the 
electrode concentrations and concentration overpotentials. We find 
that both of these quantities are proportional to a single integral as 
before, so that the electrolyte can be characterized by specifying the 
two constants of proportionality, instead of n 2 diffusivities and n 
mobilities. 
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IV. THE COMPLETE BOUNDARY VALUE PROBLEM 

The steady flow of ions in a nonreacting, concentrated, n-component 
electrolyte is governed by the n conservation laws* 

V-N,- = 0, (1) 

for % = 1, • • •, n, where the ion flux N, is given in terms of the ion 
concentrations C\, ■ ■ • , C n and the potential $ by 

N,- = dV - £ DijVCj - FziUidV*. (2) 

3=1 

In these n relations, V is the circulation velocity of the electrolyte due 
to stirring, F is the Faraday constant, z, the charge number of the 
t'th ion, and the diffusivities D,, and mobilities £/,• are functions of 
Ci, • ■ •, C n . f We shall not impose any restrictions on the form of the 
Di/s, although it seems reasonable to require, for example, that for 
i 9* j, Da — » 0, as C, — » 0. This makes the ith. ion flux N, vanish, 
when the ith. ion is absent, and gives the correct form for a dilute 
electrolyte. With 

V-V = 0, (3) 

as is appropriate for a liquid electrolyte, eqs. (1) and (2) yield the set 
of n nonlinear, second-order, partial differential equations 

V- ( £ DijVCj + FziUidV*) = V-VCi. (4) 

The set of equations is completed by the charge neutrality condition 

t Zid = 0, (5) 

t=i 

which implies that the current density 

J = E F*& 4 = -KV* - tt FziDijVCj, (6) 

»-i t-i y-i 

with the conductivity 

K = ± FhWiCi, (7) 

satisfies 

V-J = 0. (8) 



* Capital letters denote dimensional quantities ; lower case letters denote dimen- 
sionless quantities. 

* In Appendix C, we show how Z>, ; and t/, can be expressed in terms of the elec- 
trochemical potential. 
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We assume that there is no flux through the surface bounding the 
electrolyte, except that flux due to the production of cation 1 (e.g., 
Cu"*" 1 ") at the anode and deposition of the same cation at the cathode. 
Thus, for outward normal n, 

Vn = N t -n = (9) 

Ni-n - J-n/Fsi (10) 



f or i = 2, • • • , n, while 
on the boundary, where 



{—J a , on the anode, 
J c , on the cathode, (11) 

0, otherwise. 

We shall assume that the electrode current densities determined by 
the surface reaction rates have the form 

J. = J:(CyC>[exp(a e F|<f> - *.\/RT) 

-exp(-/5^|*-*.|/flr)], (12) 

for e = a, c, where the ion exchange current density J* e (Ci/C* e ) y ' is 
assumed to be a function only of the concentration of the cation 
participating in the surface reactions, R is the gas constant, T the 
absolute temperature, <i> the electrolyte potential, $, the electrode 
potential, and |$ — $«| the surface overpotential. These conditions 
complete the boundary value problem. 

Now let L, V, C, RT/F, D, D/RT, and FDC/L be some typical 
length, velocity, concentration, potential, diffusivity, mobility, and 
current density and set 

(x, y, z) = (X, Y, Z)/L, 
d = d/C, 
if, = F&/RT, 
v = V/V, 
j = LJ/FDC, 
da = Dii/D, 
Ui = RTUi/D. 

The equations and boundary conditions then take the dimensionless 
form 

V Y £ dijVcj + ZiUiCtfA = Pe v Vci, (13) 

for i = 1, • • •, n, where the Peclet number Pe = VL/D, 
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£ ZiCi = 0, 



v-(kv<p+ E E»A/Vcil = o, 



where 



k = Y, z\uid. 



With vn = on the boundary, 

n 

5Z dijdcj/dn -f- ZiUiddip/dn = 

3-1 

where 



— yn/zi, for i = 1, 
0, for i = 2, 



(14) 
(15) 

(16) 

(17) 



— j-n = kd<p/dn + £ £ ZidijdCj/dn 
t=i >=i 

on the anode, 

on the cathode, 

otherwise, 

and, for e = a, c, 

j« - itc7*[exp(o; e | 95 - p.|) - exp(-0 e | ^ - V«|)], 
/ e = (LJ:/FDC)(C/C:)y. 




(18) 

(19) 
(20) 



V. THE BOUNDARY LAYER APPROXIMATION 

To obtain explicit results, we consider a plane, rectangular region 
with a cathode at Y — y = and an anode &t Y = L, y = 1. We 
assume that the electrolyte is well stirred, with Pe » 1, so that the 
concentrations are nearly constant, except in thin boundary layers, 
of thickness of order cL, e = (Pe) _i , balancing diffusion and convection 
near the electrodes. Far from the electrodes, we introduce the outer 
expansions 



d = c { {x, y, «) = ci 0) + 0(e 3 ), 

<p = <p{x, y, e) = <p m (x, y) + e& (1) (x, y) + ■■■, 

since eq. (13) implies that v-Vc,- = 0(e 3 ), where the Ci 0> 's are the con- 
stant bulk concentrations and <p m , £ (1) are harmonic functions; i.e., 

-(0) I -^) -(1) I — tW (\ 

fix + tPm ~ <Pxx -T ipyv = u - 
In the boundary layers, we set 

d = Ziix, y, e) = cj 0) + ecP(x, y) + • • •, 
tf> = $(x, p, e) = <p w (x, y) + ■*«(«, P) + 



(21) 



(22) 
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where the stretched variable y = y/e in the cathode boundary layer 
and y = (1 — y)/e in the anode boundary layer. The assumption that 
the electrode current densities are bounded has already been intro- 
duced by choosing the leading term in the expansion of Zi to be con- 
stant. In the boundary layer eqs. (13) and (15) become 

I L ctiiCjv + ZiUiCi(f)v\_ — (ife* — $xCiy) 

= -*(£> due,* + ZiHiCi? x ) , (23) 

fei+EI Zidacji) = -<?(k<p x + t t zdijc jx ) , (24) 

\ t-iy-i II \ »=i i=i /x 

where we have introduced the reduced stream function 

$(x, y, e) = e-hl/(x, eg), 
with the x and y components of v given by 

u = *v, v = -**, 

vanishing on the electrodes. The tildes on da, etc., mean a>a 
= dij(ci, • ■ •, c„), etc. Note that the d's, ti's, k, and $ must all be 
expanded around e = to obtain the equations satisfied by ci l) , £ (0) , 
<p a) , •••. The boundary conditions, for example at the cathode 
y = ■g = 0, have the form 

V< 3 7 J- „*7 -- I *3o/Zl, for i = 1, /nrN 

Z difri + MAft =| 0) for i = 2,-,n, (25) 

where 

3c = j*e(ci) Te [exp(a fl 0) - exp(-/3 fl <p)], (26) 

Z\ and <p > being evaluated at y = and the cathode assumed to be 
grounded, so that <p e = 0. Obviously j e must also be expanded in 
powers of e. Finally, multiplying eq. (25) by Zi and summing yields 

n n 

&££ + £ £ ZidijZjy = ej c . (27) 

f-i j-i 

The set of relations from which the boundary layer equations are ob- 
tained is completed by matching Ci(x, g, e) with c,(x, y, e) and <p(x, #, e) 
with <p(x, y, e) for small y and large #. Specifically, we use intermediate 
matching, setting 

y = € "y* ) £ = e*-y, 

for fixed ?/* and < a < 1, and considering the limit e — > 0. Matching 
?,(*, e-V, «) = c t (0) + ec^(x, «-V) 
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with 

Si(x, e°y*,e) = C ( (0) + O(e 3 ) 

trivially yields the boundary conditions 

tP(fr • ) = 0, 

for i = 1, • • •, n. Similarly, matching <£ (0) and £ (0) gives 

<P w (x, oo ) = $w(x,0). 

The differential equation (24) and the boundary condition (27) 
imply that <ftj = fyf* = 0, so that 

#«(**) = <p«»(x,0). 

To a first approximation, the outer potential £ (0) , as yet undetermined, 
simply penetrates the boundary layer. 

VI. FIRST-ORDER MATCHING: THE SECONDARY CURRENT DISTRIBUTION 

Boundary conditions on £ (0) are found by matching 

with 

v = #<« + e0«>. 

The zero-order inner potential <p (0) has already been determined. The 
first-order potential £ (1) satisfies the differential equation 

L t = iy=i Ji/ 

and the boundary condition 

»-i i-i 

at # - 0, where d#>, fc< 0) are evaluated at c{ 0) , • • •, c< 0) and j™ at ci°\ 
(o, (x, 0) = £ <0) (x,0), so that 

ji 0) (s) = i![cP>{exp[a^»>(x,0)] - exp[-,W°>(x,0)]}. 

Two integrations yield 

#«(«, ji) = 0<i>(x) + [gti 0) (x) - £ .£ zdifc^(x, «]/*», (28) 

where $ (1) (x) is an arbitrary function of x. The last term can be 
identified as the concentration overpotential. Since, for large fj, c\ l) — * 0, 

$ ~ v m (x, 0) + ■[#«(*) + e-VJJ»(*)/*«l 
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for # = t a ~ l y* and e — » 0. On the other hand, expanding in Taylor 
series gives 

<p ~ 0<°>(ar, 0) + t°y*v™(x, 0) + e^Oc, 0) + • • •, 

for ?/ = e a y*. First-order matching then yields the zero-order boundary 
condition 

*"(», 0) = J2»(z)/k«» 
and 

0<»(aO = #«(*,<>), 

another penetration condition. A similar analysis produces the 
boundary condition 

at the anode y = 1. Since # 0) and ji 0) depend only on £ (0) , we now 
have a complete boundary value problem for <p w . To a first approxi- 
mation then, the potential <p satisfies Laplace's equation 

<Pxx + <pw = 0, (29) 

in the electrolyte, and the nonlinear boundary condition 

{j a , on the anode, 
- j c , on the cathode, (30) 

0, otherwise, 

where, for e = a, c, 

je = j:[c{ 0> ]^[exp(a e | <p - <p.\) - exp(-j9.| <p - *>.|)], (31) 

yielding the so-called "secondary current distribution." 

Although we have derived this secondary current distribution 
specifically for small boundary layer thickness and moderate electrode 
current density, it probably is valid over a much wider range. The 
basic assumption is simply that the concentration variations are small. 
Thus, it is probably still valid when the flow is turbulent and the 
present simple description of the flow field does not suffice or at loci 
of boundary singularities (sharp edges, etc.) where the gradients are 
not necessarily normal to the boundary, as is assumed tacitly when 
one introduces the stretched variable y in the boundary layer approxi- 
mation. In any event, by calculating the first-order terms c 4 (1) , £ (1) , 
we can obtain an error estimate or an improved approximation. 

VII. CONCENTRATION VARIATIONS 

Although we have already derived boundary conditions on the c^'s, 
we have not yet displayed the equations which they satisfy. With 
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constant leading terms c t (0) and $ expanded in the form 

= #„ + 4i + • • • = y 2 f(x) + eg 3 g(x) + • • •, 

where f(x) = ip yv (x, 0)/2, g(x) = f vvv (x, 0)/6, • • •, eq. (23) becomes, 
to first order, 

t Wt% + w /**) fS = ^ - *&*> < 32 ) 

for « = 1, • • •, n, where the ith partial conductivity 

fcf°> = z 2 M 0) cP, 
while the boundary condition (25) at y = yields 

v* j(ov~w i /mow \,<« \je m /zh for t = 1, , Q q\ 

E<^ + <WV«)* =| o, for z = 2, -..,n. (33) 

We can use the charge neutrality condition 



and the relation 



z? = - E (2,A»)^ 






found from eq. (28), to eliminate c a) and £ (1) . We obtain the reduced 
boundary layer equations 

Z4fi$l-iFl&-f&$, (34) 

for i = 1, • • •, n — 1, and the boundary conditions at # = 

V w < -w „ 1 (1 - tfOiJV*. for i = 1, .g-v 

A diiCj ' 5 " 1 -tf 8 *.»/*, for » = 2, • • •, n - 1, w 

where # 0) = ki m /k m is the ith transference number and the reduced 
diffusivity matrix [d' t j] has elements 

d' ti = dff - (ziMdW - ti 0) t M*dld& ~ (*//*-)*! ( 36 ) 



p=l 



for i = 1, • • • , n — I, j = 1, • • • , n — 1. The boundary value problem 
is completed by the conditions at infinity 

ZP»(z, oo) = 0, (37) 
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and the "initial condition" 



cP((U) =0. 



(38) 



As usual, the boundary layer approximation changes the equations 
from elliptic to parabolic type, requiring a condition at x = in the 
time-like variable x. Equation (38) is equivalent to the assumption 
that the concentrations have their free-stream, bulk values near the 
leading edge x = y = 0. We now have a complete set of coupled 
equations and boundary conditions for the first-order boundary layer 
concentration variations cl l) . We note in passing that the corresponding 
set for the second-order variations c t (2) would be much more com- 
plicated, containing not only $ (1) , as in the case of the dilute binary 
electrolyte, but also first-order terms 3$, etc., due to the concentration 
dependent diffusivities and mobilities. 

VIII. THE VECTOR CD-EQUATION AND ITS SOLUTION 

Equation (34) and boundary condition (35) can be written in the 
vector form 

Zd'J? (l) lVV = *f P a) ]- - & 0) [Z (l) ]s 
in g > 0, and 



(39) 
(40) 



OT = 



on g = 0, where [d'~] is the (w — 1) X (n — 1) matrix 

dn • • • Ol(n-l) 

"(n-l)l ' * ' «(n-l)(n-l) 

and [c (1) ] and [T 7 ^] the column vectors 
HP 



B?«] = 



r (1) , 

^n — 1 



[!T< >] = 



(1 - h)/zi 

— tn-l/Zn-l 



By way of comparison, for the binary electrolyte (n = 2), one has 

zg = fftp - m™, (4i) 

for $ > 0, and 

cf = TMj?\ (42) 

for# = 0, where ct n = 1,? 0) = VF, 7 7(0) = T[ 0) (1 - h)/z v The solution 
of this problem is given in Appendix A, together with that of the 
corresponding axially symmetric problem, important for jet plating 
and for the rotating disk electrode. 
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We now introduce the eigenvalues d {p) and the eigenvectors [w (p) ] 
of [d'], such that* 

Id'Jju^l = d (p) Ctt (p) ], 

for p = 1, • • •, n — 1. If we set 

p-i 
[T<«] = £ t< p >|><*>], 

eqs. (39) and (40) are satisfied if 

<*<*><# = *f «*■ - Wc?, (43) 

in > 0, and 

dWf = r<^j; o) , (44) 

on # = 0, for p = 1, • • •, n — 1. 

Now, in Appendix A we show that, if c(x, y) satisfies 

Cyy = fyC x ~ \{/ x Cy, 

with \p = y 2 f(x), and 

c(0, y) = c(x, oo) = 0, 

then 

c(x, 0) = [' c v (x', 0)K(x, x')dx', (45) 

Jo 

where the kernel 

K(x, x') = - (£)» (/J [/(2)]^)~y r(|). (46) 

Thus, 

c"»(x, 0) = [d^-]-*T^ [' j™(x')K(x, x')dx', (47) 

with a factor (d (p) )*, found by replacing / by //d (p) in K, and a factor 
(d (p) ) _1 from the boundary conditions. Finally, to first order in e, 
the concentrations at the cathode are given by 

Ci(x, 0) = cf°> + mt [' jFWKcfr, x')dx', (48) 

Jo 

for i ' = 1, • • •, n — 1, where 

m = E [d^H-V*^, (49) 



* We assume that such u's and d's exist. This must be verified in each particular 
case. 
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and we now write K e for K to indicate that the kernel depends on the 
velocity distribution near the cathode. A similar expression 

~ Ci (x, 0) - c t (0) + *, J* j™(x')K a (x, x')dx', (50) 

with x measured from the anode leading edge and i = 1, ■ ■ • , n — 1, is 
valid for the concentrations at the anode. The concentration of the 
?ith ion at cathode and anode is found from the charge neutrality 
relation 



n-l 



c n (x,0) = - E (ZiMti(x,0). (51) 

With the Ci (1) 's determined, we are almost ready to calculate the first- 
order corrections to the electrode current densities. For example, at 
the cathode we have 

jc = j* e (ci) yc lexp(a c <p) - exp(-/3£)], 

where 

c, = c{°> + <ci l \ 
(p = *>«» + e<p«\ 

all quantities being evaluated at y = 0. With 

j e = ti 0) +ej? ) + •••, 

we find 

ii 0) = j:[ C { 0) >[exp(a c £<°>) - exp(-/3 fl e (0) )], 

as before, and 

J c (1> = ilc! 01 ]^^/! 11 + a c <p^exp(« e <p™) 

- M 1J -p e <pVlexp(-(3 <p™)}, 

where 



n n 



#0) = ea>( Xj 0) - Z L z&MK*> o)A (0) - (52) 

We still must determine the first-order outer potential <p {l) . It is the 
solution of the linear potential problem 

9S + pf? = 0, (53) 

in the electrolyte, 

j^, on the anode, 

k w d& a) /dn = - - J e (1) , on the cathode, (54) 

0, otherwise, 

where the J c (1) 's are linear functions of <p a) , with coefficients which 
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depend on c[°\ $ m , and cl l \ This boundary condition is found from 
second-order matching of q> and $, just as the boundary condition 
on £ (0) was found from first-order matching. As is usual in boundary 
layer problems, the complexity of matching increases rapidly with 
order so this rather messy matching is relegated to Appendix B. Once 
£< l) has been determined, the electrode current densities can be calcu- 
lated to first order in e. 

IX. THE CONCENTRATION OVERPOTENTIALS 

The concentration variations cP appear in the J e (1), s in two ways. 
First of all, there is the power law dependence of the ion exchange 
densities on the concentration Z\ of the "active" cation (e.g., CU++). 
This concentration is given by 

Zx = c[ 0) + e M i P j!°KJx? (55) 

Jo 

in terms of the bulk concentration c[ 0) and the constant in. 

The second concentration dependence is contained in the so-called 
"concentration overpotential," the difference between the potential 
<p{x, 0) at the electrode and the potential $(x, 0) at the edge of the 
boundary layer. To first order in e, this concentration overpotential 
is given by 

A<Pe = e (0(D - £U>) = ed J* ji 0) K e dx', (56) 



where 



n n— 1 



5 = - L E zldif - ( Z y/2n)d 4 ( n O) >i//c (0) . (57) 

t=l }-l 



X. FORMAL PERTURBATION SOLUTION 

Although we shall eventually abandon this approach in favor of an 
equivalent scheme more convenient for numerical analysis, let us 
briefly summarize the successive steps by which the electrode current 
densities are calculated to first order in e. There are three such steps : 

(i) The nonlinear potential problem. First, we calculate the 
harmonic function <p m , satisfying nonlinear, mixed-boundary 
conditions. As is usual in boundary layer theory, we begin with 
a nonlinear problem, the only nonlinear problem in the se- 
quence. It yields first approximations to the electrode current 
densities, which may be adequate for many purposes. It also 
furnishes a boundary condition for the boundary layer con- 
centration variations. 
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(ii) Concentration variations. The concentration variations c< (1) , or 
rather ?j 1} and A<p e , which are the only quantities involved in 
J e (1> , are calculated from eqs. (55) and (56), giving them in 
terms of a single integral of the product of electrode current 
density and given kernels. 
(Hi) The linear potential problem. With c[ l) , A<p« known, the har- 
monic function <p a) , satisfying linear, mixed-boundary condi- 
tions in £ (1) , c[ l) , A ^ e is calculated and, from it, a second ap- 
proximation to the electrode current densities, good to order e. 

XI. SOLUTION BY BOUNDARY INTEGRALS 

Instead of following the above rather cumbersome procedure, we 
can adopt the approach suggested by J. L. Blue, 2 forming a nonlinear, 
coupled set of integral equations whose solution, e.g., by successive 
approximation, yields the electrode current densities directly to order 
«. From this point of view, the formal perturbation scheme simply 
validates this method. For example, perturbation tells us that the 
integral solution of the boundary layer equation is good to order e 
and no more. 

We begin with Green's third identity for a plane region : 

2*<p(P) = f [j a (Q)/fc]G(P, Q)ds(Q) - f [i.(Q)/*] 

J T a JTe 

■G(P, Q)ds(Q) - J^ V (Q) |? (P, Q)d(Q), 

where G has the fundamental singularity at P = Q, T a is the anode, 
r„ the cathode, and r the whole boundary, and we now denote £ (0) 
+ e£ (1) by <p and the bulk conductivity k m by k. For each set of 
values of j a , j B , (p, and points P on the boundary, the solution of this 
integral equation gives <p. Similarly, from the integrals 



C, = c{°> + m J* j e K e dx', 

A<p e = +e8 [ j e K e dx', 
Jo 



with e = a, c, we can calculate c x and A<p e from the given values of 
j a and j c . With <p, ci, and A^„ now known on the electrodes, we can 
calculate new values of j a and j e from the relations 

je = j*e(Cl)""[exp(ae\ <p + A ?« - tp,\) ~ exp(-/3 e | <p + Aipe ~ <P«|)], 

for e = a, c, the anode potential <p a > 0, and the cathode potential 
(p c = (for a grounded cathode) . Using these new values of electrode 
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current densities, compared with the old values, a successive approxi- 
mation scheme can be developed which converges rapidly. 

Note that this scheme is formally identical with that for the dilute, 
binary electrolyte. Furthermore, within the framework of this calcu- 
lation, the concentrated, n-component electrolyte can be characterized 
completely by giving its bulk conductivity k and the constants n\ and 
8. The values of the individual diffusivities da and mobilities w,-, from 
which jui and 8 can be calculated, are of no importance. Only k, m, 
and 8 enter the final set of equations. 

APPENDIX A 

Integral Solution of the Boundary Layer Equation 

Let c(x, y) satisfy the boundary layer equation 

Cm = tvCx — fxCy (58) 

in x > 0, y > 0, with ^ = y 2 f(x) and 

c(0,y) = c(x, co) =0. (59) 

For use in the method of boundary integrals, we wish to express the 
boundary values of c in terms of an integral over the boundary values 
of c y , i.e., we seek the kernel function K in the relation 

c(x, 0) = [* c y (x', 0)K(x. x')dx'. (60) 

Jo 

Following Levich 3 and Lighthill, 4 we make the change of variable 

s = f* imydz, v = gw» = ytmy, 

the form of rj being motivated by the observation that the right-hand 
side of eq. (58) vanishes when c is a function of $ alone. We find that 
V (Z, v) = c(x, y) satisfies the equation 

»« = 2ifv f , (61) 

if £ > 0, 77 > 0, and 

v(0, v) -»(«, ») =0. (62) 

Denote the boundary values of v v by g, so that 

»,(*,<>) = g(Q. 
Then the Laplace transform v, given by 



Hv, s) = J" e~* y(£, ij)dt, 



BOUNDARY LAYER 819 



satisfies the ordinary differential equation 

v" - 2sr,v = 0, 

in < v < oo , and the boundary conditions 

S'(0,s) =g(s), «J(oo, s) =0, 

where g is the transform of g. Now suppose w(£, i?) is any solution of 
eq. (61), satisfying eq. (62). Then we have the identity between 
transforms 

Hv, *) = g(s)u(v, s)/u'(Q, s), 

giving 

v(0, s) = f»'(0, «)[tf(0, s)/w'(0, a)]. (63) 

A convenient choice for u is the similar solution of eqs. (61) and (62), 
satisfying the condition 

«(*, o) = i, 

so that w(0, s) = 1/s. In this case, u is a function of f = *?/£* alone, 
given by 

M = [/ t "exp(-2 2 V9)<fe]/(i)'r(i), 

so that 

«,(€, o) = -rV(£)*r(£) 

and, with T(k)/s k the transform of £ 1_fc , 

«'(o,s) = -r(|)/(i)*r(i)A 

Thus, 

«(o, «)/«'(o, «) = -(i)*r(|)/r(f)s*. 

Inversion of eq. (63), using the convolution theorem, then yields 
where k, the inverse of w(0, s)/m'(0, s), is given by 

«(0 = - (l)*r f /r(f). 

Using the relations 

»(€', 0)d£' = c„(a;', 0)dx', 

€ - r = f c/w: 1 *, 

./ a' 

we finally obtain eq. (60), with 

K M = -(!)'(/; C/«]^)7r(|). (64) 
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K is of convolution type, i.e., a function of x — x' alone, only when 
/ is constant, e.g., for Poiseuille flow near the wall of a channel or 
tube. For the laminar boundary layer flow on a semi-infinite flat plate, 

/ ~ ar*, 

so that 

K ~ (x* - x'*)-». 

A slight modification of the above yields the integral relation 

c(r,0)= f ' c,(r', 0)L(r, r')r'dr', (65) 

Jo 

for the axially symmetric problem, where 

C tz = T-^ypzCr - frCz), 

in r > 0, z > 0, with ^ = 2 2 /(r) and 

c(0, z) = c(r, co) = 0, 
where the kernel function 

*fr.o --(£)*(£ WH^VKi) - (66) 

We merely replace the previously defined £ by 

./o 

Equations (65) and (66) apply, for example, to the rotating disk elec- 
trode and to "jet plating," where a jet of electrolyte is directed against 
an electrode. 

APPENDIX B 

Second-Order Potential Matching 

To obtain a boundary condition for the first-order outer potential 
£ (1) , we must match 

!p = £(°> + ev m + «V 2) 

with 

q> = £ (0) + e£ (1) + eV 2) , 

where 

p(o) = ^(O) fo 0), (67) 

#« = ^a)( Xj o) + *;•/*« - EL *.<W/* (0) - (68) 
The second-order inner potential <£ (2) satisfies the equation 
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and the boundary condition at # = 

*»# + £«>fc a) + ii 1.(41^; + <wS) = jp, (70) 

where k, da, j c have been expanded in the forms 

k = /c< 0) + ek^ + •••, 
da = iff + «Zg> + •••, 

ic = ir + «jP»+ •••• 

Note that £ (1) , (Iff, j c (1) are linear functions of the ci^'s; for example, 
Two integrations give 

0O> = 0W + J|J<»/*<W - ^<°>/2 _ £ f 0i ^(O)^(2) 

t-1 j-1 



+ 



r [ 1 1 **tf ° w/*»> - WHS - k^jw/kw i ^/*», 

./o l»-i y-i J 

where $ (2) is an arbitrary function of x and eq. (68) has been used to 
eliminate <p a) . We wish to calculate <p (2) in the limit # = e 0-1 ?/*, with 
< a < 1, y* fixed, e— > 0. The first three terms all contribute to 
the limit, while the fourth term vanishes, since c (2) — > 0, as $ — > °o . 
Finally, we must evaluate the limit of the integral. Consider the last 
term first. Since $ 0) is a function of x only and 

|d) - £ (dk«»/dc P )c?\ 

p-i 

we have n integrals like 

Jo 
Changing variable gives 

f" c^(x, $)dg = e*" 1 f c™(x, e a - l z)dz-*e^ 1 I(e), 

Jo Jo 

where 1(e) — > 0, as c— »0. Similar limits are obtained for the other 
integrals with integrands of the form c^c^. Thus, 

~ £ (0 >(z, 0) + O a) («, 0) + e-ViT/JfeW] + e 2 [^«>(x) 

+ e a-i y * <7 -(D //c (0) _ ( e «-y*) 2 ^°>(z, 0)/2], 

since pfi? = $2£(x, 0), to be matched with 
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<5~ *»>(«, 0) + e«y*ip?\x, 0) + (e a y*) 2 <P$!(x, 0)/2 

+ «[>«(*, 0) + .V^fo 0)] + **«(*, o). 

Since £ (0) is harmonic, $£?(£, 0) = - $S(x, 0), and matching term by 
term gives the previous boundary condition on £ (0) 

<p?\x, o) = j«v*<°>, 

the boundary condition on £ (1) 

and the "penetration" condition 

0»>(a;) = ^<W( ajj 0). 

APPENDIX C 

The Electrochemical Potential 

In Section IV, we expressed the fluxes in terms of the gradients in 
the form 

N,- = CiV - t DiFCj - FzifrUiV*, (71) 

in effect defining the diffusivities D 1; and mobilities U i} functions of 
the concentrations, by this relation. In electrochemistry, it is more 
usual to give the gradients of the electrochemical potentials, i.e., well- 
defined thermodynamic functions, in terms of the fluxes. Thus, New- 
man 1 gives the relations (for n + 1 species in our notation) 

ftVw = t KiWt - V,), (72) 

for i = 0, 1, • • •, n, where the m's are the electrochemical potentials, 
the Ki/s are the interaction coefficients, functions of the concentra- 
tions, and V.'s are subject to certain constraints. Thus, 



»=o 
and 



L Win = o 



£ p.v,- = pV, 



t=0 



where V is the fluid velocity, p,- = Mid is the partial density for 
molecular weight M if and 



P = H Pi- 

t'-0 
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These restrictions imply that 

± (K t , - K 3i ) = 0, 

»'=0 

for j = 0, 1, ■••, n. In particular, Newman assumes the Stefan- 
Maxwell relations Ka = Ka, i.e., pairwise interaction between species. 
With these restrictions, eq. (72) has a unique solution for V,- = V. 
Thus we find, following Newman, that 

V t - V = - E L ip C p Vn P , 

for i = 1, • • • , n, where — [Lif] is the inverse of [_M ,-,•], with 

n 

- E #»fc, for * = j, 

Ka, for i ?^ j. 

(Note that we may assume that Ka = 0, for i = j.) Once V,- — Vo is 
calculated, V* — V may be determined from the identities 

Vo - V = - E (p*/p)(V* - Vo), 

k=i 

V, - V = V, - Vo + Vo - V. 
With 



V/*p = E (dn P /dCj)vCj + (<Wd*)v$, 



we finally find 

Do, = - II (i> k /p)L kl fW,(dn,/dC,), 

A«=l p-1 

^of/o =-II (pML kp c p (d^ p /e^) t 

A=l p=l 

fl« - E |~£<p - E (p*/p)£*,1 dC^dpJdCi), 

p = l L * = 1 J 

ft«l7< = E [l,p - E (pML kp ] c p (<Wa*). 

P =i L *=i J 
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